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Abstract. We argue the possibility of realization of a polarization-insensitive all-optical switching in a 
planar metamaterial composed of a 4-fold periodic array of two concentric metal rings placed on a substrate 
of nonlinear material. It is demonstrated that a switching may be achieved between essentially different 
values of transmission near the resonant frequency of the high-quality-factor Fano-shape trapped-mode 
excitation. 
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■ 1 Introduction the bistable device is a Fabry- Perot interferometer, filled 
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with a Kerr nonlinear material. In this case, the resonator 
Optical bistability is an area of modern optics which is provides feedback, which is essential to obtain a multival- 

•• ; evolving rapidly. The great attention to this phenomenon ue d intensity at the structure's output. However, in such 

^ , 
• i— i 

^ ; is due to the possibility of its practical application to de- a system, both relatively strong light power and/or large 
H ■ 

03 sign optical switches, limiters, transistors and diodes. Typ- enough volume of nonlinear optical material are generally 
ically, in bistable devices, the nonlinear medium is placed needed to achieve a sizeable nonlinear response, 
inside an optical cavity, just as is done in lasers, but unlike 
the latter they are passive bistable devices whose operat- 



To overcome these drawbacks, photonic crystal micro- 
cavities [5] and quantum well structures [3] were proposed 

ing conditions are easy to control because the light prop- . ,. ,. „ „ . 

to enhance the nonlinear crtccts as well as to reduce the 



agation is controlled with light [J . A classical example of 



material volume. With the assistance of surface plasmon 
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cal optical field intensity, optical bistability has also been mode excitation of an actual structure, the held is strongly 

shown numerically in different metal nanostructures such localized to the structure plane and resonant transmission 

as surface plasmon polaritonic crystals [1] , subwavelength and reflection have a large quality factor due to very small 

gratings that consisted of infinitely long slits in metal radiation of electromagnetic energy in a comparison with 

slab (5], jB], etc. In all these cases, the excitation of high- stored one. 

quality- factor resonances in the systems is provided to ob- • ln , • i u - , , j i 

lypically, metamaterials wmcn can bear trapped modes 

tain efficient switching. r-j , i i_ i ^ ±. n- • i 

° consist of identical subwavelengtn metallic inclusions struc- 

A promising way to produce optical switching in com- tured in the form of asymmetrically split rings [12] , |13) , 
pact devices can be found also in using planar metamate- split squares |14) , |15j or their complementary pattern |16) , 
rials. It is known that the planar metamaterials can create |15) . These elements are arranged periodically and placed 
an environment equivalent to a resonant cavity. As usual on a thin dielectric substrate. In the structure of such kind, 
such structures are composed of metallic elements in the the high-quality-factor current oscillations with the lowest 
form of symmetrical split-ring resonators [7] , [8] . They are total emission losses appear when all currents in the metal- 
resonant because of an internal capacitance and indue- lie elements of a periodic cell oscillate in antiphase. The 
tance within each element. Most often such structures are characteristic feature of such metamaterials is the depen- 
investigated in order to obtain a negative refraction index dence of their spectra on the polarization and the angle of 
over a finite frequency range, although their nonlinear op- incidence of input waves. In |17j . |18j . the polarization in- 
tical response was also studied [S]-[TT]. Unfortunately a sensitive structure configurations also were proposed. One 
considerable disadvantage of such structures in the con- such structure consists of a planar array with a periodic 
text of obtaining optical switching consists in their low cell element that consists of two concentric rings (double- 
quality factor. ring (DR) structure). 

Nevertheless, exceptionally strong and narrow reso- In the DR-structure, the trapped-modes excitation is 

nances are possible in planar metamaterials via engag- exclusively controlled by the difference in circumferences 

ing trapped modes }12) . The quality factor of the system of the inner and outer rings and does not depend on the 

depends on ratio of power of stored energy to power of ra- polarization state of the normally incident electromagnetic 

diation and dissipation losses. If we suppose infinitesimal radiation. In contrast to the asymmetrically split rings 

dissipation losses in the structure, the trapped modes cor- metamaterials, DR-structure is a system, where mutual 

respond to real eigenvalues (i.e., real resonant frequencies) coupling between elements of neighboring cells is weak 

of the relevant boundary value problem due to a special ge- and the response of the entire array is practically a direct 

ometry of metal elements. In the regime of quasi-trapped sum of the DR individual contributions. As a consequence, 
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the electromagnetic response of the DR-metamaterial is 
weakly dependent on the angle of incidence |19] . Remark- 
ably, at the trapped-mode resonance the electromagnetic 
energy is confined to a very small region between the rings, 
where the energy density reaches substantially high val- 
ues. This makes the response of the metamaterial oper- 
ating in the trapped-mode regime extremely sensitive to 
the dielectric properties of the substrate. This feature can 
be used for enhancing optical nonlinear response in the 
nanoscaled version of the metamaterial [30] ■ 

The goal of this paper is to show promising use of pla- 
nar metamaterials that bear the trapped-mode resonances 
to obtain all-optical switching. Especially, we argue the 
possibility of realization of a polarization-insensitive all- 
optical switching in planar metamaterial designed on the 
basis of double-ring array placed on a substrate of nonlin- 
ear dielectric in the regime of trapped-mode excitation. 

2 Problem statement and solution 



Fig. 1. (Color online) Fragment of the planar metamate- 
rial and its elementary unit cell. 

index of these materials in the mid-IR range is between 
1CT 4 cm 2 /kW and 10" 2 cm 2 /kW. 

As the excitation field a normally incident plane monochro- 
matic wave of a frequency lu and an amplitude A is se- 
lected. We suppose that the intensity of the incident field 
is enough for the nonlinearity to become apparent, i.e., it 
is about 1 kW/cm 2 . 

In the set of our previous works [H], [T3], [H]-[2D] the 
method of moments |25| was used to calculate the response 



The square unit cell (d — d x = d y ) of the structure of different planar metamaterials in the microwave range 

under study consists of one DR (Fig. [TJ. The radii of when the amplitude of the incident wave is small (linear 

the outer and inner rings are fixed at oi/d = 0.36 and case). This numerical method involves solving an integral 

a 2 /d = 0.29, respectively. The width of both the metal equation for the surface currents induced in the metallic 

rings is 2w/d = 0.05. The array is placed on a nonlin- pattern by the incident electromagnetic wave, then calcu- 

ear dielectric substrate with thickness h/d = 0.2. As a lating the scattered fields produced by the currents as a 

suitable variant of substrate material we propose semicon- superposition of partial spatial waves. In the framework 

ductors appropriate in the mid-IR region of wavelength, of this method it is implied that the metallic pattern is a 

for example, InSb or InSb with some impurities such as thin perfect conductor. 

As, Tl, Bi, P [21]. From the experimental results [22]- [24], On the other hand in [15], on the basis of the pseudo- 



we expect that the Kerr-nonlinear part of the refractive spectral time-domain algorithm, the optical response of 
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the metamaterial in the IR range was calculated taking as 

into account the strong dissipation and dispersion of the Ei n — V/b, 

metal permittivity of the elements. It was revealed that where y = ZJ ig the Une vo i tagei b = ai _ fl2 _ 2w , 

the trapped-mode resonances are well observed and have j is the ma g nitude f current wh i c h flows along the DR- 

high Q-factor in this range, and the obtained results are in element) and z is the inl p ec i ance f i ine . The impedance is 

good agreements with the ones evaluated with the method determined at the resonant dimensionless frequency ae = 

of moments down to the mid-IR region. rf/Ao 

Thus in the present paper we will use the method of g _ ggif^2_ 

dCo 

moments to calculate the magnitude of the current J along 



the single ring of the DR element, the reflection r and 



where I = n(ai + o.2)/2, and 
1 



C o = -ln 

transmission t coefficients. They can be determined in the 4 



JL + JfJLY-! 

2w V \2wJ 



f orm: is the capacity in free space per unit length of line, p = 

ai — a-i- From this model it follows that the electric field 



J = J(uj,e), t = t(uj,e), r = r(uj,e). (1) 



Next we suppose that the structure substrate is a Kerr 
nonlinear dielectric whose permittivity e depends on the 
intensity of electric field i, n inside it. Under rigorous con- 
sideration, the nonlinear permittivity e of the substrate 
is inhomogeneous. The permittivity riches its maximum 
value directly under the metallic pattern and along the 
rings this permittivity is also different. Nevertheless, as 
mentioned above, at the trapped-mode resonance, the elec- 
tromagnetic energy is confined to a very small region be- 
tween the rings and the crucial influence of the permittiv- 
ity on the system properties occurs in this place. There- 
fore, the approximation based on the transmission line The incident field magnitude A is a parameter of equa- 



strcngth between the rings is directly proportional to the 
current magnitude J. Since the unit cell is small in com- 
parison with the wavelength, the current magnitude J can 
be substituted with its value averaged along the metallic 
ring, J. 

Thus, from this model it follows that the electric field 
strength between the rings is directly proportional to the 
average current magnitude J, and the nonlinear equation 
on the average current magnitude in the metallic pattern 
is obtained in the form 



J = A-J(LJ,£ 1 +e 2 Ii„(J)). (2) 



theory can be used here to estimate the field intensity tion (0). At a fixed frequency w, the solution of this equa- 

between the rings. According to this theory, conductive tion is the average current value which is dependent on 

rings are considered as two wires with a distance b be- the magnitude of the incident field, J = J (A). 
tween them. Along these wires the currents flow in oppo- On the basis of the current J(A) found by a numerical 

site directions. Thus the electric field strength is defined solution of equation @, the value of the permittivity of 
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the nonlinear substrate e = e± + E2li n {A) is determined 
and the reflection and transmission coefficients (H} are cal- 
culated 

t = t(u,si+ e 2 Iin{A)), r = r(w,£i + £ 2 I in (A)), (3) 

as the functions of the frequency and magnitude of the 
incident field. 

It should be noted here that our treatment of nonlin- 
earity in the planar metamaterial differs from the ones 
considered earlier in [5]-[TT] where some effective medium 
parameters were introduced to obtain the nonlinear dou- 
ble negative materials. 

3 Numerical results 

In the case of the nonlinear permittivity of substrate, de- 
pendences of the inner intensity and the transmission co- 
efficient magnitude versus the incident field magnitude 
lin = lin(A) are typical and have the form of hystere- 
sis (Fig. [5]). Such form of curves of the transmission coeffi- 
cient magnitude is studied quite well [T] and can be under- 
stood from the following considerations. Suppose that the 
trapped mode resonant frequency is slightly higher than 
the incident field frequency. As the intensity of the inci- 
dent field rises, the magnitude of currents on the metal ele- 
ments increases. This leads to increasing the field strength 
inside the substrate and its permittivity as well. As a re- 
sult, the frequency of the resonant mode decreases and 
shifts toward the frequency of incident wave, which, in 
turn, enhances further the coupling between the current 
modes and the inner field intensity in the nonlinear sub- 
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Fig. 2. (Color online) (a) The inner intensity and (b) 
the magnitude of the transmission coefficient versus the 
incident field magnitude in the case of the nonlinear per- 
mittivity of the substrate. For this and further calcula- 
tions we use the substrate parameters E\ = 4.1 + 0.02i 
and £2 = 5 x I0 -3 cm 2 /kW. The value of dimensionless 
frequency ae = d/X is chosen a bit lower to the frequency 
of the trapped-mode resonance (aeo = 0.3). 

strate. This positive feedback increases the slope of the 
rising edge of the transmission spectrum, as compared to 
the linear case. As the frequency extends beyond the reso- 
nant mode frequency, the inner field magnitude in the sub- 
strate decreases and the permittivity goes back toward its 
linear level, and this negative feedback keeps the resonant 
frequency close to the incident field frequency. As a result, 
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Fig. 3. (Color online) Frequency dependences of (a) the 
inner intensity and (b) the magnitude of the transmission 
coefficient in the case of the nonlinear permittivity (ei = 
4.1 + 0.02i and e 2 =5x 10~ 3 cm 2 /kW) of the substrate. 



ities. The Fano-type resonance appears as a result of the 
interference between a high-quality resonance and a much 
smoother, continuum-like spectrum and typically exhibits 
a sharp asymmetric line shape with the transmission coef- 
ficients varying from to 1 over a very narrow frequency 
range |26j . Such form of resonance is very suitable to ob- 
tain great amplitude of switching since there are gently 
sloping bands of the high reflection and transmission be- 
fore and after the resonant frequency. 

4 Conclusion 

In conclusion, a planar DR nonlinear metamaterial, which 
bears a high-quality-factor Fano-shape trapped-mode res- 
onance, is promising object for a realization of a polarization- 
insensitive all-optical switching. 

This work was supported by the National Academy of 
Sciences of Ukraine under the Program " Nanotechnologies 
and Nanomaterials" , Project No. 1.1.3.17. 



at a certain intensity of the incident field, the transmission 
coefficient stepwise changes its value from small to large 
level. 

The frequency dependences of the transmission coef- 
ficient magnitude also manifests discontinuous switching 
from small to large level with frequency increasing (Fig. |3j) . 
This switching appears closely to the resonant frequency 
of the trapped-mode excitation. The main peculiarity of 
the observed resonance is that this trapped-mode reso- 
nance in DR-structure has the Fano-shape rather than the 
Lorentzian one, as is the characteristic of Fabry-Perot cav- 
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